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Abstract 

We consider calibration of the Black-Karasinski short rate model to a given interest rate (or credit 
intensity) term structure, conducting an asymptotic analysis in the limit of low rates. We base our 
analysis on the perturbation expansion approach proposed by Turfus and Schubert (2016). We calibrate 
the model so as to give consistent representation of the probability distribution of T maturity zero coupon 
bond prices at all times up to maturity. 


1 Introduction 

As pointed out recently by Turfus (2016b), the Black-Karasinski lognormal volatility short rate model is 
known to have many attractive features as an interest rate model; however its most significant perceived 
drawback since it was first proposed 25 years ago has been its lack of tractible analytic formulae for conditional 
zero coupon bond prices or caplets. This was recently addressed by Turfus (2016b) using a perturbation 
expansion based on an assumption of small volatility. It was proposed from the form of the expansion 
derived that, even when the lognormal volatility is not particularly “small”, the caplet formulae ought still 
to be reliable, provided that the rates themselves are not too large. 

It is often however wished to use the Black-Karasinski interest rate model as a credit intensity model, 
given its ability to ensure positive rates. This feature makes it more attractive than the otherwise popular 
(because relatively tractible) normal volatility short rate model of Hull and White (1990). But it is the 
square root model of Cox, Ingersoll and Ross (1985) which tends to be most favoured for credit intensity 
modelling because it is both tractible and can be configured to ensure positive rates. None the less, as has 
been pointed out by Brigo and Mercurio (2006), it too has a significant downside as a credit model insofar 
as it is limited in its ability to calibrate to markets with implied lognormal volatility in excess of about 
40-50%. For this reason we propose below a perturbation expansion approach to allow calibration of a 
Black-Karasinski model to a given forward curve (interest rate or credit intensity) subject to the assumption 
that the level of the forward curve remains low, but without any assumption about the level of the volatility. 

Having set out our model and modelling assumptions in §2, we proceed by calculating in §3 the theoretical 
distribution under our model of T-maturity zero coupon bonds at arbitrary future times t. By applying a 
condition of the consistency of the implied present values with the market-specified forward curve, we are 
able to tie down the model to a unique calibration in an asymptotic sense of errors in bond prices being 
formally O(e)!), where e r is a parameter capturing the “smallness” of the short rates under consideration. 1 It 
is proposed that, even for fairly stressed credit intensity rates, this ought to provide sufficient accuracy for 
many practical purposes, in particular risk management, which is where a great deal of the focus is in credit 
modelling in the finance industry at present. 

1 Or, equivalently, the relative error in the calculated rates = O(e)l) 
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2 Modelling Assumptions 


We consider short rate processes representing the interest short rate to be driven by a Black-Karasinski 
model. We shall find it convenient to work with a reduced variable y t satisfying the following canonical 
Ornstein-Uhlenbeck process: 

dy t = - a r y t dt + a r (t) dW t (1) 


where dW t is a standard Brownian motion under the martingale measure. This is taken to be related to the 
interest short rate r t by 

r t = r*(t)e yt- 5 / >'( t o. t ) ( 2 ) 


where ^ 

I r (ti,t 2 ) := [ e~ 2o ‘ r( ' t2 ~ s ^a 2 (s) ds. 

Jt i 


(3) 


We take to to be the initial observation date with, by assumption, y to = 0. The function r*(t) is determined 
by calibration but is constrained for consistency to tend in the zero volatility limit to the forward interest 
rate r(t). The formal no-arbitrage constraint which determines this function is as follows 




r(s) ds 


(4) 


under the martingale measure (money market numeraire) for to < t < T m , where T m is the longest maturity 
date for which the model is calibrated. We proceed by writing the (stochastic) price at time t of a unit 
cash flow at time T as We then look to price these cash flows and so to determine the conditions 

necessary to satisfy Eq. 4. 

For future notational convenience, we also introduce at this stage the deterministic discount factors 

D(t 1 ,t 2 )=e-fn r ^ ds (5) 


3 Calibration of Model to Zero Coupon Bonds 


We deduce by standard means that the V maturity zero coupon bond price fr{y,t) will be governed under 
the money market numeraire by the following backward diffusion equation: 


dfr 

dt 


- ot r y 


dfr 

dy 


+ Writ) 


d 2 fr 
dy 2 


r*{t)e v f T , 


( 6 ) 


the r.h.s. being simply fa multiplied by our assumed functional representation of r t . 

In the absence of exact closed form solutions to Eq. 6, we seek an approximate solution under a low 
interest rate assumption, i.e. we rescale f{t) by a non-dimensionalised asymptotic parameter e r defined by 


fl7 r W dt 

^r(^ra t 0 ) 


and a scaled forward rate 

f(t) := e“ 1 r(t) 

where f(t) is 0(1) as e r — >• 0. On this basis we expect the following asymptotic form for r*(t): 

r*(t) = e r r{(t) + e 2 r r* 2 {t) + 0(e 3 r ) 


(7) 


We can then rewrite Eq. 6 as 


£[My,t)\ = (e r (rt(t) 


e »-pr(to,t) _ f ( t ) 


+ elr. 


r'2\ 


(t)ey-^*))f T (y,t) + 0(4) 


( 8 ) 


2 



where £[.] is a standard forced diffusion operator given by 2 

£[1 = I “ + “ F<<) ’ 


and the final condition that must be satisfied is 

Mv,T) = 1. 

To solve for fr(y,t) we pose an asymptotic expansion 

fr(y,t ) = /o(t,T) +e r /i((/,t,T) + elf 2 {y 7 t,T) + 0(e 3 r ). 
At zeroth order we have that £[/o(i)] = 0 with the trivial solution that 

fo(t,T) = D(t,T), 

which clearly satisfies the required final condition. At first order we have 


£[fi(y,t,T)] = - f(t)) f 0 (t,T) 


The final condition Eq. 10 requires that = 0. The no-arbitrage condition Eq. 4 introduces the 

further constraint that fi(0,to,T) = 0. 

To proceed we observe that Eq. 9 has a Green’s function solution given by 

G(y,t;r],v) = D(t,v)H(v - t)^-N (- — ^ \ (13) 

dy \ y/I r {t,v) J 

where H{.) is the Heaviside step function and N(.) is a standard unit normal distribution. We deduce by 
standard means: 

fi{y,t,T) = -D{t,T) j exp (e~ ariv ~ t] y - \e^ 2o ‘ r{v ~ t) I r {t 0l t))J - r(vf) dv. 

Applying the no-arbitrage condition fi(0,to,T) = 0 we conclude we must choose 

rt(t) = r(t) (14) 

whence our first order result can be written in terms of the Doleans-Dade exponential function £(.) as 


fi(Vt,t,T) = -D(t,T)J T (£ (e- aAv - t] yt) - l) r(v), 


Proceeding in a similar vein, we find at second order 


C[f 2 (y,t,T)}=f(t) ( e »-* J -(W) - l) My^ + rWW-^’Vfofrn (16) 

with solution given by 

h{yt,t,T ) 

rT rex) 

= -J J G(y t ,t-,y,v) r(v) (e r, -^ Ir ( to ’ v) - lj /i(r?,?i,T) + r2{v)e n ~i I ^ to ’ v) f 0 (v,T) dydv 

= D(t,T) J T f(u) £ (£ (e-^-Vyt) £ (e~ aAv -^y t ') 6 xp (e"“- ^I r (t,vj) 

- £ (e~ otr ^ v ~ t ) y t ^j - £ + ^ dv du 

- D(t,T) J T r* 2 (u)£ (e-“-(«-t)y t ) du 

2 Note that we have kept r(t) in its original unsealed form within the operator definition since this drives only the discount 
factor scaling, in which context we will not seek to expand in terms of the small parameter. 
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This clearly satisfies the final condition f 2 (yr,T,T) = 0. At t = to we have 


f 2 (0,t 0 ,T) = D(t 0 ,T) J (r(u) ^ F*(v,u)r(v)dv - r 2 (u)^ du 


where we have defined 

F*(v,u) := exp (^e~ ar( - u ~ v ^ I r (to, v — 1 
The initial condition / 2 (C), icb^) = 0 then requires us to choose 


( 17 ) 


(18) 


r 2 (u) = r(u) f F*(v,u)r(v) dv. 

Jt 0 

Finally we can write 3 

h(Vt,t,T) = \D{t,T) ^ (£ (e-^-Vyt) - l) f(t;)d^ 

/> T p U 

+ D(t,T) J 8 [e~°‘ T ^ u ~ t ' , yt S j f{u) J F*(v, u) ^8 (^e~ ar ^ v ~ t ^yt^ — 1^ r(v) dv du 
— D(t,T) J 8 [e~ aT< ' u ~ t ' , yt^ r («) F*(v,u)r(v)dvdu (19) 

Combining all contributions and reverting to unsealed notation, we obtain our main result: 

Mvt,t) = D(t,T) 1 - £ (s [eT^^y^j - l) r(v)dv + \ ^ (s [eT 01 ^^ y^ - l) r(v) dv^j 
+ ^ 8 (e~°‘ T ^ u - t ^yt^r(u) J F*(v, u) [8 (e~°‘ r ^ v ~ t ^y t ^ — 1^) r(v) dv du 
— 8^e~°‘ r( ' u ^ t:> yt S jf(u)J^ F*(v,u)r(v) dv du 


o(4) 

The other main results are Eqs. 14 and 18 which in unsealed notation can be written 

r*(u) = r(u) 

r 2 (u)=f(u) f F*(v,u)r(v) dv. 

J t 0 

This gives the model calibration, accurate to O(e^). 


( 20 ) 
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